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We report electron spin resonance (ESR) measurements in the Gd 3+ doped semiconducting filled 
skutterudite compound Cei_ a ;Gd I Fe4Pi2 (x ~ 0.001). As the temperature T varies from T ~ 150 
K to T ~ 165 K, the Gd 3+ ESR fine and hyperfine structures coalesce into a broad inhomogeneous 
single resonance. At T ~ 200 K the line narrows and as T increases further, the resonance becomes 
homogeneous with a thermal broadening of 1.1(2) Oe/K. These results suggest that the origin of 
these features may be associated to a subtle interdependence of thermally activated mechanisms 
that combine: i) an increase with T of the density of activated conduction-carriers across the In- 
dependent semiconducting pseudogap; ii) the Gd 3+ Korringa relaxation process due to an exchange 
interaction, Jf^S.s, between the Gd 3+ localized magnetic moments and the thermally activated 
conduction-carriers and; Hi) a relatively weak confining potential of the rare-earth ions inside the 
oversized (Fe2Ps)4 cage, which allows the rare-earths to become rattler Einstein oscillators above 
T w 148 K. We argue that the rattling of the Gd 3+ ions, via a motional narrowing mechanism, also 
contributes to the coalescence of the ESR fine and hyperfine structure. 



I. INTRODUCTION 

The filled skutterudite RT4X 12 compounds, where R is 
a rare-earth or actinide, T is a transition metal (Fe, Ru, 
Os) and X is a pnictogen (P, As, Sb) have attracted great 
attention due their broad range of physical properties. In 
particular, they are of interest to those investigating basic 
mechanisms of strongly correlated electronic systems^^ 
and also to those seeking for more efficient thermoelectric 
materials!^ 

These compounds crystallize in the LaFe4Pi2 structure 
with space group ImZ and local point symmetry T/ t for 
the R ions. The R ions are guests in the oversized rigid 
^2X3)4 cagesPThe dynamics of the guest R ions is be- 
lieved to be of great importance in the damping of the 
thermal co ndu ctivity observed in the filled skutterudite 
compounds LSISI Moreover, they may also play an impor- 
tant role in the appe aran ce of heavy fermion behavior 
and superconductivity!^ 

Electron spin resonance (ESR) is a sensitive and pow- 
erful microscopic tool that provides information about 
crystal field (CF) effects, site symmetries, valencies of 
paramagnetic ions, <?-values, and fine and hyperfine 
parameters P^l I n a recent work our groups found ESR 
to be a sensitive and useful tool to study the dynamics 
of the R ions in this family of filled skutterudites. The 
weak confining potential on the R ions at the center of 
the oversized cage allows them to easily get off-center 
and experience slightly different local strength and sym- 
metry of the CF which may lead to: i) a distribution of 
the ESR parameters and, ii) a rattling of the R ions that, 
due to motional narrowing effects, 12 may cause remark- 



able changes in the observed ESR spectra. In our pre- 
vious ESR experiments on Yb 3+ in Cei_ a Yb a: Fe4Pi2p^ 
these two features were observed and the coexistence of 
two distinct Yb 3+ sites was confirmed. 

Ogita et aZpSl performing Raman scattering experi- 
ments on several metallic skutterudite compounds of the 
RT4X12 (T = Fe,Ru,Os; X = P,Sb) series, found reso- 
nant 2 nd order phonon modes associated with the vibra- 
tion that change the bond length of the R-X stretching 
mode. However, in semiconducting CeFe4Pi2 the 2 nd or- 
der phonon modes were found to be non-resonant. Based 
on their results Ogita et aP^l concluded that there should 
be a strong coupling between the R-X stretching modes 
and the conduction-electrons (c-e). Most reports on the 
T-dependence of the rfc-resistivity in CeFe4Pi 2 present 
a semiconductor- like behavior!^ However, the resistiv- 
ity is strongly sample dependent and only in some cases 
it shows metallic behavior below T « 200 K.--- : Never- 
theless, for most of the reported samples the conduc- 
tivity due to thermally activated carriers predominates 
above T rj 200 K. Thus, for the semiconductor CeFe4Pi2 
with gap of ~ 0.15 eV and estimated Debye tempera- 
ture of 0£> ~ 500 K, 46 at least a weak coupling of the 
R-X stretching mode and the c-e should be expected. 
This compound experiences a huge increase in the den- 
sity of thermally activated conduction-carriers at T w 
150 K. 14 Also, evidences for rattling of the Yb 3+ and 
Ce 4+ ions were found in ESRP^ and extended x-ray ab- 
sorption fine structure (EXAFS) 1 ^ experiments, respec- 
tively. The aim of this work is to learn if the presence 
of thermally activated conduction-carriers and rattling of 
the R ions can be observed by the ESR technique. For 
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that reason we measured the evolution of the Gd 3+ ESR 
spectra in Cei_ a: Gd 2 ;Fe4Pi2 with T. To compare our 
data with a nonrattling compound, we have also studied 
the evolution of the Gd 3+ ESR spectra in Cai-^Gd^Bg 
(x ~ 0.001) with T, which is a cubic CsCl type semicon- 
ductor with a gap of ~ 0.8 eV and Debye temperature of 
e D ~ 783 Kpl 

We found that for Cei_ x Gd :E Fe4Pi2 the ESR spectra 
show a different behavior in three T-regions. At low- 
T the system behaves as an insulator, at high-T as a 
metal and in the intermediate region it presents the ef- 
fects of: a) an exchange interaction, J/^S.s, between 
Gd 3+ localized magnetic moments and thermally acti- 
vated conduction-carriers and 6) possible evidence for 
rattling of the R ions. 



II. EXPERIMENTAL 

Single crystals of Cei_ a ,Gd a; Fe4Pi2 (x < 0.001) were 
grown in Sn-flux as described in Ref. HU The cubic 
structure (7m3) and phase purity were checked by x-ray 
powder diffraction. Crystals of ^2x2x2 mm 3 of natu- 
rally grown crystallographic faces were used in the ESR 
experiments. Single crystals of Cai^Gd^Bg (x < 0.001) 
were grown as described in Ref. [TH1 The cubic structure 
(space group 221, Pm3m, CsCl type and local point sym- 
metry Td for the R ions) and phase purity were checked 
by x-ray powder diffraction and the crystals orientation 
determined by Laue x-ray diffraction. Most of the ESR 
experiments were done in ~ 2x1x0.5 mm 3 single crys- 
tals. The ESR spectra were taken in a Bruker X-band 
(9.48 GHz) spectrometer using appropriated resonators 
coupled to a T-controller of a helium gas flux system for 
4-2 < T < 300 K. The Gd concentrations were deter- 
mined from the H and T-dependence of the magnetiza- 
tion, M(H,T), measured in a Quantum Design SQUID 
dc-magnetometer. In both systems the magnetic suscep- 
tibility follows a Curie-Weiss behavior. Also, in both 
compounds the T-dependence of the Gd 3+ ESR intensity 
presents a Curie- Weiss like behavior within the accuracy 
of the experiments. 



III. RESULTS AND DISCUSSION 

In both compounds at low-T the Gd 3+ ESR spectra 
show the full resolved fine structure corresponding to the 
spin Hamiltonian for the Zeeman and cubic CF inter- 
actions, H — g/3HS + 64O4 + b^Oe™ The angular and 
T-dependence of the spectra were taken mostly with the 
applied magnetic field, H, in the (1,-1,0) plane. The 
fitting of the data to the spin Hamiltonian shows that 
the parameters are, within the experimental accuracy, 
T-independent for the entire studied T-range. The mea- 
sured parameters were, g = 1.986(3) and 64 = 7(1) Oe 
for Ce 1 _ x Gd a; Fe4Pi 2 and g = 1.992(3) and 64 = 13.8(5) 
Oe for Cai-^Gd^Bg, in agreement with previous low-T 
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Figure 1: T-dependence of the X-band ESR spectra for H 
[001]; a) Cei-.Gd^P^ and b) Cai-^Gd.Be. 



reports 1221211 The accuracy of the data was not enough 
to estimate the value of 6 6 . The g-shift measured for 
Gd 3 + in Cei_ x Gd a; Fe4Pi2 is negative, Ag = 1.986(3) - 
1.993 « - 0.007. An additional term to the spin Hamil- 
tonian, J/rfS.s, due to a covalent exchange hybridization 
between the Gd 3+ 4/ electrons and c-e with d charac- 
ter would be responsible for this negative g-shift. 22 For 
the Cei_ a; Gd ;E Fe4Pi2 crystal, careful measurements of 
the spectra were taken from T ~ 150 K to T ~ 200 
K for various directions of H. In this T- interval, the 
fine structure coalesces into a single broad line and its 
lineshapc changes from lorentzian (insulator) to dysonian 
(metallic) Notice that these features are independent 
of the field orientation and none of them is observed in 
Cai-^GdzBg. 

Figures 1, 2 and 3 display the evolution with T (4.2 
< T < 300 K) of the normalized ESR spectra of Gd 3+ 
in the Cei_ x Gd a ,Fe4Pi 2 and Ca^^Gd^Bg (x < 0.001) 
crystals for H in the (1, — 1, 0) plane along [001], 30° from 
[001], and [110], respectively. These data show that for 
T > 150 K the T-dependence of the Gd 3+ fine structure 
is quite different in both compounds. 

For Ce 1 _ x Gd a; Fe4P 1 2 the central transition (| «-» — i) 
at 4.2 K and H along [001] is narrow enough to observe 
the hyperfine satellites lines of the isotopes 155 > 157 Gd 3+ 
(I = 3/2) (see Fig. 4). The measured hyperfine parameter 
is A = 5.5(2) OePS This hyperfine structure is also ob- 
served, although not so clearly, for the other transitions 



3 



Ce^Gd/e/,, X = 0.001 9.5 GHz H // - [30' 



A- 



Ca, Gd B x = 0.001 



-v 

-V 



3.3 3.4 3.5 

H (kOe) 



Figure 2: T-dependence of the X-band ESR spectra for 
H along 6 « 30° from [001] in the (1,-1,0) plane; a) 
Cei_ a: Gd ;c Fe4Pi2 and b) Cai-^Gd^Be. 
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Figure 3: T-dependence of the X-band ESR spectra for H 
[110]; a) Cei_,Gd a; Fe4Pi2 and b) Cai-^Gd.Be. 



in the spectrum. For the angle where the fine structure 
collapses (29.6° from [001]) and the various transitions 
overlap, a small missorientation of H by < 2° away from 
this direction affects the overall lineshape and the hy- 
perfine structure is then strongly blurred (see Fig. 2). 
For Ca^^Gd-rBg Fig. 4 shows that, due to its higher 
(/-value and broader linewidth, the (^ o — ^) transition 
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Figure 4: X-band ESR (| O - §) transition for H || [001]. The 
arrows show the satellite hyperfine structure for both samples 
of x = 0.001. 



is shifted to lower- H and the hyperfine structure is not 
well resolved. However, the hyperfine parameter can still 
be estimated to be A ~ 7(1) Oe. 

Figure 5 presents for both compounds the evolution 
with T (4.2 < T < 300 K) of the linewidth (AH) for the 
various Gd 3+ ESR transitions at several H orientations. 
For Cei_ a; Gd 2 ,Fe4Pi2 Fig. 5a shows the T-dependence of 
AH for the (| o -i) transition and H along [001], [111], 
[110] directions and at - 30° from [001] in the (1, -1,0) 
plane for the collapsed spectrum. It is clear from the data 
that there are three regions of different T-dependence 
of AH: Region I, for T < 150 K, where AH is nearly 
T-independent and very narrow at ~ 5(1) Oe; Region 
II, for 165 < T < 200 K, where the full fine structure 
dramatically coalesces into a broad inhomogeneous single 
resonance with anisotropic AH: and Region III, for T > 
200 K, where AH is again isotropic and homogeneous. 
It corresponds to a single coalesced resonance and has a 
linear thermal broadening of ~ 1.1(2) Oe/K, reminiscent 
of a Korringa-like relaxation process via the c-eP^ 

Figs. 5b and 5c show the T-dependence of AH for the 
various transitions in Cei_. E Gd. E Fe4Pi2 and Cai-^Gd^Bg 
for H along [001] and 30° from [001], respectively. A 
timid broadening starts to be observed on AH for the fine 
structure components at T ~ 60 K for Cei_ 2; Gd a ;Fe4Pi2 
and at T ~ 120 K for Cai_ x Gda;B6. Presumably this 
broadening is caused by a phonon spin-lattice relax- 
ation processPS The fact that such a phonon contri- 
bution starts at lower-T in Cei__ x Gd :E Fe4Pi2 than in 
Cai-zGd^Bg is consistent with the lower Debye temper- 
ature for the former compound. Alternatively, the Gd 3+ 
ions produce bound states in the gap, which in the case 
of CaBg are donor states. Carriers bounded at low-T in 
these states can be promoted into the conduction band 
as T increases and produce a faster relaxation. How- 
ever, as T increases in Region I a small local distribution 
of CF cannot be excluded as the reason for the small 
broadening of the fine structure lines. The large voided 
space and concomitant increase of the carrier density as 
T increases may thermally activate slow motions of the 
Gd 3+ ions inside the oversized (Fe2Ps)4 cage which could 
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Figure 5: T-evolution (4.2 < T < 300 K) of the Gd 3+ ESR 
linewidth, AH, for both compounds and various transitions 
at several H orientations. Notice that in Fig. 5a the behavior 
of AH(T) clearly characterizes three different Regions (I, II, 
and III). The solid lines in Regions II and III correspond the 
calculated AH(T) for the coalesced ESR spectra using the 
Plefka-BarneP^I exchange narrowing mechanism (see text). 
Fig. 5d presents the general T-dependence reported for the 
resistivity in these compounds. The resistivity of our crystals 
is similar to that of Ref. 14. 



slightly alter, in an inhomogeneous way, the local CF at 
the Gd 3+ site. 

As already mentioned, for Cei_ :r Gd. E Fe4Pi2 above T ~ 
160 K a dramatic broadening mechanism drives the whole 
Gd 3+ resolved ESR fine structure in Region I to coalesce 
into the broad inhomogeneous and unresolved anisotropic 
spectrum of Region II (see Fig. 5a). This striking result 
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the field for resonance between the low and high-T spectra is 
mainly due to the change, from lorenztian to dysonian, and 
also to a small change in the frequency of the microwave cavity 
due to the temperature (<0.05%). 



occurs at about the same T where: i) the density of ther- 
mally activated mobile carriers increases by several orders 
of magnitude (see Fig. 5d); it) the rattling of the filler 
R atom is confirmed by EXAFS experiments^ Hi) the 
existence of a T-dependent semiconducting pseudogap is 
observed for T < 300 K in ultraviolet and x-ray photoe- 
mission spectroscopies (UPS, XPS)p3and; iv) where the 
change from lorentzian (insulator) to dysonian (metallic) 
ESR lineshape is observed (see Figs, la, 2a and 6). Note 
that none of these features are present in Ca^^Gd^Bg. 
The solid lines in Regions II and III of Fig. 5a are the 
calculated AH(T) for the coalescing ESR spectra using 
the Plefka-Barnes^ exchange (J/^S.s) narrowing theory 
of the fine structure. In the calculation we used a Kor- 
ringa relaxation of 1.1(2) Oe/K that is ' 'switched- on" at 
157(2) K, a fourth order CF parameter 64 = 7(1) Oe and 
a residual linewidth AH(T = 0) = 5(1) Oe. 

Figure 6 presents the T-dependence of the hyperfine 
structure for the O —5) transition. The data show 
that the coalescence of the hyperfine structure is already 
observed at T w 150 K, i.e., at ps 15 K below the coa- 
lescence of the fine structure at T rj 165 K. This is ex- 
pected since the exchange interaction would act first on 
the hyperfine structure due to its much smaller spectral 
splitting. 

Figure 7 shows the angular dependence of AH at dif- 
ferent T corresponding to Regions II and III. Following 
the analysis of Urban et aZpSl for the exchange narrowing 
of the Gd 3+ ESR fine structure, the anisotropy of AH in 
Region II can be fitted to the general expression for the 
intermediate coupling regime: 



where 



AH = A(T) + B(T)p 2 (9), (1) 



p 2 (0) = 1 - 5 [sin 2 (<9) + (3 /4) sin 4 (0)] . (2) 
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Fig. 5a shows that there is narrowing of AH for T ap- 
proaching 200 K and that the anisotropy decreases, i.e., 
B -> as T -> 200. For T > 200 K Ai? becomes isotropic 
(see Fig. 7) and increases linearly at a rate of 1.1(2) 
Oe/K (see Fig. 5a). This linear increase is an evidence 
for the presence of a Korringa relaxation process, i.e., the 
Gd 3+ ions relax to the lattice via an exchange interac- 
tion, JfdS.s, between the Gd 3+ localized magnetic mo- 
ment and the thermally activated conduction-carriersP^ 

A <?-shift of Ag = —0.007 has been measured for the 
entire range of T studied for Cei_ I Gd x Fe4Pi2. This is 
surprising because in Region I there are no conduction- 
carriers that could be polarized. However, the host is a 
Kondo insulator with a finite Van Vleck susceptibility due 
to the crystalline field splitting of the Ce ions. This Van 
Vleck susceptibility is larger than the susceptibility of the 
thermally excited electrons in Region II and III and pro- 
vides the polarization to produce the (/-shift. This effect 
is of course not present in the Cai^Gd^Bg sample, since 
CaB 6 has no significant susceptibility. There is a second 
unusual issue with the g-shift. In a simple metallic host 
Gd 3+ ions are expected to have a ferromagnetic Heisen- 
berg exchange. However, the g-shift is negative, indica- 
tive of a hybridization mechanism. The overlap of the 
Gd 4/ electrons with the hybridized Ce 4/-band forming 
the valence and conduction bands of the Kondo insulator 
could give rise to an antiferromagnetic exchange. 

It is interesting to note that in Gd 3+ doped simple met- 
als the Korringa relaxation, d(AH)/dT, would be related 
to the g-shift, Ag, byPl 
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Figure 7: Angular dependence of the AH for the Gd 3+ ESR 
for temperatures in the Regions II and III of Fig. 5a. The 
solid line corresponds to the fitting of the data to Eq. 1 for 
A(172 K) = 10(3) Oe and B(172 K) = 80(10) Oe. 
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Using our experimental value of 1.1(2) Oe/K for 
d(AH)/dT and 2.34 x 10 4 Oe/K for nk/gfi B we esti- 
mate a corresponding Ag of ss |0.007|. These results 
and Eq. ^ suggest that: i) in Region I, where there 
is no Korringa relaxation, the exchange coupling due to 
covalent hybridization gives rise to just polarization ef- 
fects, J/d(q = 0)j^3 ii) the trigger of the Korringa mech- 
anism in Regions II and III is due to the presence of 
mobile activated conduction-carriers at the Fermi level, 
which are responsible for the momentum transfer be- 
tween the conduction-carriers and the localized magn etic 
moment via the exchange coupling, J/d(q 7^ 0)j 3 ° l 31 l and 
in) in the metallic Regions II and III, there is no q- 
dependence of the exchange interaction, i.e., J/d(q = 
0) EE <j| d (q))V™ 

For T < 10 K Figure 8 shows that, due to the long 
spin-lattice relaxation time T\ in these materials, the 
Gd 3+ (| o — i) transition saturates as a function of 
the microwave power P^l From AH = ^X^) -1 the spin- 
spin relaxation time can be estimated to be T 2 ~ 0.1/xs 
in both compounds. From the saturation factor, S = 
[l+(l/4)i/ 1 2 7 2 TiT 2 ]- 1 and P = (l/4)iff, we estimate 
the spin-lattice relaxation time to be T\ ~ 10 ms and 
T\ ~ 4 ms for Cei_. E Gd. E Fe4P 12 and Cai-^Gd^Bg, re- 
spectively. Notice that a Korringa relaxation is absent 



Figure 8: Dependence of the Gd' j+ ESR intensity at low-T on 
the microwave power for both materials. 



in Region I where the compound behaves as an insula- 
tor. Moreover, from AH at T ~ 300 K in Fig. 5a we 
estimate T\ ~ 0.002/is which is much shorter than the 
low-T value of T 2 . Therefore, at high-T T 2 ~ T x and, 
as far as ESR is concerned, this is another evidence that 
Cei_ x Gd :E Fe4Pi 2 behaves as a regular metal in Regions 
II and III. 

Finally, according to the Raman results and 
conclusions about the rattling modes in metallic 
skutteruditesj^ it is plausible that the huge increase in 
the metallic character of CeFe4Pi 2 activates the R-X 
stretching mode and triggers the rattling of the R ions 
inside the oversized (Fe 2 P3)4 cage. Hence, via a mo- 
tional narrowing mechanisnP^ the rattling of the Gd 3+ 
ions could also contribute to the dramatic change of the 
ESR spectra observed at the transition from Region I to 
Region II. Actually Figure 9 shows that the exchange nar- 
rowing mechanism alone cannot reproduce the observed 
experimental single coalesced resonance at the transition 
between these two Regions. Thus, a motional narrow- 
ing of the Gd 3+ fine structure, due to rattling of the R 
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Figure 9: Experimental and calculated 111 Gd 3+ ESR spectra 
in Cei_ a; Gd a: Fe4Pi2 for x = 0.001 at the transition from Re- 
gion I to Region II. The parameters used to calculate the 
theoretical spectra were the same as those used for the calcu- 
lation of AH(T) in Fig. 5a. 

ions, cannot be disregarded and it probably contributes 
to the experimental observed spectra. It should be men- 
tioned that such a striking behavior is not expected in 
Cai-^Gd^Bg due to the much larger semiconducting gap 
of s» 0.8 eV and the tighter cages for the CaBg compound. 

IV. CONCLUSIONS 

In this work we have presented experimental data that 
show the following striking features for the T-dependence 
of the Gd 3+ ESR spectra in Cei_ x Gd 2: Fe4Pi2: a) the 
coalescence of the hyperfine and fine structures at T ~ 
150 K and T ~ 165 K, respectively; b) at about these 
temperatures the ESR lincshape changes from lorcntzian 



(insulating media) to dysonian (metallic media); and c) 
the T-dependence of the ESR AH changes from a narrow 
nearly T-independent linewidth for each fine structure in 
Region I to a single inhomogeneous broad resonance with 
anisotropic AH in Region II and then, in Region III, to a 
homogeneous linewidth with a broadening which is linear 
in T, resembling the Korringa-like relaxation process in a 
metallic hostPSPoint b) indicates that, at our microwave 
frequency and between T ~ 150 K and T ~ 165 K, there 
is also a clear and strong change in the ac-conductivity 
of the material. We associate this change to a smooth 
crossover from insulator to metal which was only possible 
to be detected due to the high sensitivity that the ESR 
lineshape has in a metallic media. 

Our ESR observations in Cei_. E Gd :c Fe4Pi2, along with 
those of Raman, EXAFS, UPS and XPS for CeFe 4 Pi 2 , 
suggest that this smooth insulator-metal crossover may 
be responsible for the coalescence and narrowing of the 
hyperfine and fine structures and also for the activation 
of the R-X stretching mode that probably triggers the 
rattling of the Gd 3+ ions in the oversized (Fe2Pa)4 cage. 
Via a motional narrowing mechanism the Gd 3+ rattling 
may also contribute to the dramatic change of the ESR 
spectra at the transition from Region I to Region II j 1 1 1 12 l 

We believe that our ESR study gives further clues and 
insights for the subtle interplay between the local vibra- 
tion modes (Einstein oscillators) of the R ions and the 
c-e in the filled skutterudite compounds. In particular, 
our work supports the idea that some metallic character 
is always needed to set up the necessary conditions for 
the rattling of the R ions in these materials. 
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